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Abstract 

A d-dimensional quantum model in the spherical approximation confined 
to a general geometry of the form L'^~^' xoo^' xL^ (L-linear space size and 
I/^-temporal size) and subjected to periodic boundary conditions is consid- 
ered. Because of its close relation with the quantum rotors model it can be 
regarded as an effective model for studying the low-temperature behavior of 
the quantum Heisenberg antiferromagnets. Due to the remarkable opportu- 
nity it offers for rigorous study of finite-size effects at arbitrary dimensionality 
this model may play the same role in quantum critical phenomena as the pop- 
ular Berlin-Kac spherical model in classical critical phenomena. Close to the 
zero-temperature quantum critical point, the ideas of finite-size scaling are 
utilized to the fullest extent for studying the critical behavior of the model. 
For different dimensions 1 < d < 3 and < d' < d a detailed analysis, in 
terms of the special functions of classical mathematics, for the free energy, 
the susceptibility and the equation of state is given. Particular attention is 
paid to the two-dimensional case. 
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I. INTRODUCTION 



In recent years there has been a renewed interestoQ in the theory of zero-temperature 
quantum phase transitions initiated in 1976 by Hertz's quantum dynamic renormahzation 
groupi for itinerant ferromagnets. Distinctively from temperature driven critical phenomena, 
these phase transitions occur at zero temperature as a function of some non-thermal control 
parameter (or a competition between different parameters describing the basic interaction 
of the system), and the relevant fluctuations are of quantum rather than thermal nature. 

It is well known from the theory of critical phenomena that for the temperature driven 
phase transitions quantum effects are unimportant near critical points with > 0. It 
could be expected, however, that at rather small (as compared to characteristic excitation 
in the system) temperature, the leading T dependence of all observables is specified by the 
properties of the zero-temperature critical points, which take place in quantum systems. The 
dimensional crossover rule asserts that the critical singularities of such a quantum system 
at T = with dimensionality d are formally equivalent to those of a classical system with 
dimensionality d + z {z is the dynamical critical exponent) and critical temperature > 0. 
This makes it possible to investigate low-temperature effects (considering an effective system 
with d infinite space and z finite time dimensions) in the framework of the theory of finite- 
size scaling (FSS). The idea of this theory has been applied to explore the low-temperature 
regime in quantum systems (see Refs. when the properties of the thermodynamic 

observables in the finite-temperature quantum critical region have been the main focus of 
interest. The very quantum critical region was introduced and studied first by Chakravarty 
et ali using the renormahzation group methods. The most famous model for discussing 
these properties is the quantum nonlinear 0{n) sigma model (QNLaM). 

Recently an equivalence between the QNLcrM in the limit n ^ oo and a quantum version 
of the spherical model or more precisely the "spherical quantum rotors" model (SQRM) 
was announced.0 The SQRM is an interesting model on its own. Due to the remarkable 
opportunity it offers for a rigorous study of finite-size effects at arbitrary dimensionality 
SQRM may play the same role in quantum critical phenomena as the popular Berlin-Kac 
spherical model in classical critical phenomena. The last one became a touchstone for various 
scaling hypotheses and a source of new ideas in the general theory of finite size scaling (see 
for example Refs. and references therein). Let us note that an increasing interest 

related with the spherical approximation (or large n-limit) generating tractable models in 
quantum critical phenomena has been observed in the last few years.0&ii 

In Ref. the critical exponents for the zero-temperature quantum fixed point and the 
finite-temperature classical one as a function of dimensionality was obtained. What remains 
beyond the scope of Ref. O is to study in an exact manner the scaling properties of the 
model in different regions of the phase diagram including the quantum critical region as a 
function of the dimensionality of the system. In the context of the finite-size scaling theory 
both cases: (i) The infinite (i-dimensional quantum system at low-temperatures oo'^ xL^ 

{Lr ~ ( ) finite-size in the imaginary time direction) and (ii) the finite system 

confined to the geometry L'^^'^'xoo'^'xL^ (L-linear space size) are of crucial interest. 

In this paper a detailed theory of the scaling properties of the quantum spherical model 
with nearest-neighbor interaction is presented. The plan of the paper is as follows: we start 
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with a brief review of the model and the basic equations for the free energy and the quantum 
spherical field in the case of periodic boundary conditions (Section H) . The relation of this 
model with the model due to Schneider, Stoll and BeckETEj has been briefly commented in 
Section |T|. Since we would like to exploit the ideas of the FSS theory, the bulk system in the 
low-temperature region is considered like an effective {d + 1) -dimensional classical system 
with one finite (temporal) dimension. This is done to enable contact to be made with other 
results based on the spherical type approximation e.g. in the framework of the spherical 
model and the QNLaM in the limit n — oo. The scaling forms of the free energy and of the 
spherical field equation are derived for the infinite (Section |1V A ) and finite (Section |1V B ) 
system confined to the general geometry L'^^'^ xoo'^ xL^. In Section |V| we analyze in detail 
the equation for the spherical field. This equation turns out to allow for analytic studies of 
the finite-size and low-temperature asymptotes for different d and d'. Special attention is 
laid on the two dimensional system. The remainder of the paper contains the details of the 
calculations: Appendices ^ and |B[ 



II. THE MODEL 



The model we will consider here describes a magnetic ordering due to the interaction of 
quantum spins. This has the following formlll 



H 



laT.'P'-lll -^""^A + ^J2^^-Hj2s„ (2.1) 



where are spin operators at site i, the operators are "conjugated" momenta (i.e. 
[5^,iS£,] = 0, [PgjVf^,] = 0, and [P^jS^] = iS^^,, with h = 1), the coupling constants J^^^, = J 
are between nearest neighbors only,t3 the coupling constant g is introduced so as to measure 
the strength of the quantum fluctuations (below it will be called quantum parameter), H is 
an ordering magnetic field, and finally the spherical field /i is introduced so as to ensure the 
constraint 

E ('^'> = (2.2) 

e 

Here is the total number of the quantum spins located at sites "f of a finite hypercubical 
lattice A of size Li x L2 x ■ ■ ■ x = N and (■ ■ ■) denotes the standard thermodynamic 
average taken with Ti. 

Let us note that the commutation relations for the operators and together with the 
quadratic kinetic term in the Hamiltonian ( p.l|) do not describe quantum Heisenberg-Dirac 
spins but quantum rotors as it was pointed out in Ref. |l^. 

The free energy of the model in a finite region A is given by the Legendre transformation 

fA{P,g,H) := sup In (/3, ^, i/; /i) - , (2.3) 

where Z/^{P, g, H; fi) = Tr [exp (— /5?-^)] is the partition function of the model and P is 
the inverse temperature with the Boltzmann constant Kb = 1- Under periodic boundary 
conditions applied across the finite dimensions the free energy takes the form 
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^ K q 



2 sinh ( -Puj (q; fi) 



2c<;2 (0; /i) 



(2.4) 



Here the vector q is a collective symbol, which has for Lj odd integers, the components: 

27ml 27rnd 1 \ Lj — 1 Lj — 1 



Ld J 



and 



^'^ {q; i^) = 9 yj- - '2J Y1 

The supremum of the free energy is attained at the solutions of the mean-spherical con- 
straint, Eq. (|2.2|) , that reads 



2^ 

2N 



E 



1 



+ 2E?=i(l -cosgi) 



coth 



2t 



+ 2^(1 -cosgi) + 



k2' 



(2.5a) 



1=1 



which is equivalent to the following 



t 

N 



E E 



m=—oo q 



+ 2X]i=i(l - cosgj) + b'^m'^ 



+ 



(2.5b) 



where we have introduced the notations: A = is the normalized quantum parameter, 
t = j-the normalized temperature, h = - the normalized magnetic field, b = and 
(f) = ^ — 2d is the shifted spherical field. 



Eqs. (|2.4|) and (|2.5| ) provide the basis of the study of the critical behavior of the model 
under consideration. 

A previous direct analysi^l of Eq. ( p.5| ) in the thermodynamic limit shows that there 
can be no long-range order at finite temperature for c? < 2 (in accordance with the Mermin- 
Wagner theorem). For d > 2 one can find long-range order at finite temperature up to a 
critical temperature tc(A). Here we shall consider the low-temperature region for 1 < d < 3. 
Before passing to the corresponding analysis we would like to make some comments on the 
relations of the model considered here with other models known in the literature. 



III. RELATIONS TO OTHER MODELS 

Recently&S another model, suitable to handle the joint description of classical and 
quantum fluctuations in an exact manner, was considered. If we consider Ve and Si as 
canonically conjugated momentum and coordinate, respectively, of the i atom with mass 
at each lattice point i G Z*^, then the first three terms in Eq. (|2.1| ) describe a harmonic 
interaction of the £-th and £'-th atoms. One can see that in the absence of a spherical 
constraint, if A = ^ + ^ < 0, such a lattice is unstable i.e. the parameter A < defines the 
frequency of a mode unstable in the harmonic approximation suggesting that an appropriate 
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stabilization of the lattice can create a gap in the phonon spectrum. In the spirit of the 
self-consistent phonon approximationil, it is possible to add the term 



B 
AN 



J2S^] , B>0, (3.1) 



"switching on" an anharmonic interaction which stabilizes the lattice. Because the term is 
inversely proportional to the particle number the model under consideration turns out to 
be exactly soluble in the thermodynamic limit .ilB 

In normal coordinates the "new" model Hamiltonian reads 

Knh. = IY1 (^'J^-'? + ^li<l)^9^-9) + 4^ E S,^-,^,'^-,' (3.2) 

1 qq' 

where Vg and Sg are Fourier components of the operators Vi and Si, respectively. The 
frequency of phonons is given by 



[q) = -ul + 2gJ ^ (1 - cosg,) . (3.3) 



4 2 

The anharmonicity constant h = g'^B = where Eq = is the barrier height of the 
double-well potential in (|3.2|) at a uniform displacement of all particles: Se —>■ x; 

U{x) = -^x' + ^x\ 

In Ref. |3l| the following pseudoharmonic approximating Hamiltonian has been proposed 

KppJA) = ^Y1 + ^'(9' ^)SqS-q) - ^^o(l + A)2, (3.4) 

q 

where the trial harmonic frequency 

nl{q,A)=ul{q) + vl{l + A) (3.5) 
is defined by the equation for the temperature-dependent gap A (c.f. with Eq. ( p.5|) ) : 

1 . A . -4 E ^ ^ E ^ CO.. (3.6, 

Eq. (|3.6| ) plays the role of the "soft spherical constraint". If B goes to zero, the quartic 
self-interaction disappears and our model becomes a pure harmonic model. One can see 
that this limit is singular and that 5, for d more than the upper critical dimension, will be 
a dangerous irrelevant variablcil The model (|3.2| ) is a quantum counterpart of the "soft" 
classical mean spherical model studied in Ref. ^ in the context of FSS theory. 

The inspection of Eq. ( p.6|) shows that it is similar to Eqs. ( p.5|) up to a linear term in 
the l.h.s of (|3.6|) (This term appears to be essential only above the upper critical dimension.) 
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Thus there is a clear mathematical analogy between the model defined in Section |I| and the 
model presented here. 

One of the main problems in studying the FSS effects for this type of models is to 
answer the question: how correct is the finite-size description of the initial model (|3.2| ) on 
the base of the approximating Hamiltonian ( |3.4| ). Evidently the answer to this question will 
clarify the more general and subtle problem of the status of finite-size results obtained by 
approximating methods. A successful step in this direction is presented in Ref. ^ when 
the classical mean spherical model in the Husimi-Temperely limit has been considered. In 
Ref. an appropriate modification of the approximating Hamiltonian approach, which 
reproduced some exact FSS results, has been suggested. Let us note that, if we parallel 
Ref. |3^, for studying quantum systems a more refined mathematical problems must be 
solved and that up to now this is an open problem. 



IV. SCALING FORM OF THE FREE ENERGY AT LOW TEMPERATURES 



Let us denote by /a {P,g,H) the reduced free energy /a {l3,g,H) /J . Then, from Eq. 



( p.4[ ) one obtains immediately 



/a {t, A, h) =snp{^Y^ In 



2 sinh — 
2t 



cosgj 



1=1 



(4.1) 



Using the identities 



, sinh b 1 ^ , 6^ + Tr^m^ 

^^^"1 — =o / 5-^' 

smh a 2 ^ — ^ + vr^^m^ 

m=—oo 



(4.2) 



where ab > 0, a,b are arbitrary real numbers and 



ln(a + 6)=lna+ / exp (— ax) (1 — exp (— 6a;)) — , 
Jo X 



(4.3) 



where a > 0, a + b > the above expression for the free energy can be rewritten in the form 



U{t,X,h)='- 5^ ^exp[-x(0 + 6W)] 



m=—oo 



1 - -^Sl, (x) X ... X -^SlJx) 



+tln 



2 sinh ( A 



2 6 2^ 



(4.4) 



Here 



L-l 



Sl{x) = X^exp 



m=0 



—2x 1 — cos ■ 



2Txm 

IT 



(4.5) 
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and in ( |4.4| ) is the solution of the corresponding spherical field equation (|2.5| ). In the 
remainder of this article we will be interested only in systems with geometry L'^~'^ x oo'^ x L^. 
With the help of the Jacobi identity 



oo oo / 2 \ 



(4.6) 



the above expression for the free energy can be rewritten in the following equivalent form 

d-d'' 



2/A(t,A,/i) 



A 



A 



dx 







3/2 



exp {—xq 



[exp (-2x) /o(2x)] 



d' 



dx 



X- 



3/2 



X [exp(-2x)/o(2x)]' 



exp(-x0)i?(^ 
V2t In 



d-d' 



A 



2sinh ( — 



+ 2rf)--, (4.7) 



where 



R{x) = exp (— xm^) 



(4.8) 



m=l 



and Io{x) is a modified Bessel function. 

Dividing the integrals from to cxd in the r.h.s of ( [4. 7] ) into integral from to and 
from to oo, after using for these intervals the corresponding asymptotics of ^^(x)!! 



= 1 + 2R ( -^3; ) — Lv{x) + 0(exp(— x const)) 



(4.9) 



for X > L^, where 
and 



v{x) = (Anx) ^^'^ [l — erf (ttv^)] 



Sl{x) = Lexp (-2x) /o(2x) + {L/^/^) R (^^^ 



(4.10) 



for < X < L^, one ends up with the following expression for the free energy at low 
temperatures (j ^ 1) 



2Uit,X,h) = A/rf(0) + Av^-(</. + 2rf)-- 



dx 

+A / — (47rx) 
Jo ^ 



'id+l)/2 



exp [—xq 



X 



(4.11) 
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where 



fA4>) 



I dx 
Atx Jo x^I^ 



1 - (exp(-2a;)/o(2a;))' 



The above expressions are vahd for any d. For 1 < c/ < 3 it is easy to show that 

4r(^) 



/d(0)=/d(O)+>V,(O)0 



(rf2 - 1) (47r) 



(rf+l)/2 



where 



>Vd(</.) = ^^^y c?gi---y cig, U + 2^(1 -cos 

is the well-known Watson type integral. 

Combining ( [4.11| ) and ( |4.13|) and denoting Ac := 1/Wrf(0) we obtain 



1fL{t,\h) = A/d(0)-2rf + 
+A / — (47rx 



a: 



- 1 u- 



4r (^) A 



(4.12) 



(4.13) 



(4.14) 



{d^ - 1) (47r)('^+^)/' 



L(d+l)/2 _ 



-{d+l)/2 



exp {—x(, 



1 - ( 1 + 2i? ( — 1 1 ( 1 + 2i? ^ 



\Ax J 



d~d' 



(4.15) 



A. The infinite system 

Here we consider the behavior of the free energy of the bulk system for low temperatures. 
From Eq. ( [4.15|) , taking the limit L — >■ cxd one derives immediately 



2L (t, A. k) = A/, (0) - 2. + (A _ i) ^ _ _i£glA_^(.«./. 



- 2A / — (Anx)-^^^''^' exp (-a;0) i? ( — 1 , (4.16) 



where is the solution of the corresponding equation for the bulk system (see Eq. ( ^.1|) 
below). Formally, one can obtain the corresponding equation by requiring the first derivative, 
with respect to 0, of the r.h.s of (|4.16|) to be zero. Therefore, the singular (the dependent) 



part of the free energy can be rewritten in the form 



fs^n9,bulk {t, A, h) = ^A6^+^X, {xx, , (4.17) 



where 
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1 1 



r^'^-^) (4.18) 



and 

Z^;, = /^r(^+3)/2A-i/2 (4.19) 

are the scaling variables and 

V f u \ (^) {d+i)/2 hi 

_2 jT" ^ (4»i)-<'+"'' exp (-ly,) R (^^^ (4.20) 

is the universal scaling function for the free energy. In the last expression yx := is the 
solution of the spherical field equation 



B. The finite system 



For the finite system from Eq. ( |4.15|) , in a full analogy with the previous results for the 
infinite system, we derive 



f sing, finite 

[t, A, h, L) — —Xb ~^ X finite hx, a) 



where the new scaling variable 

a = Lb 

is introduced, whereas the scaling function now is 



4r (— ) 

X finite {xx, hx, a) = -xxyx - ^^2 _ (4^)(^+l)/2 ^^' 

+ [ — (Anx)^^'^^^''^'^ exp {—xyx) 
Jo X 



(d+i)/2 _ hi 

yx 



1 - 1 + 2i? 



X 



(4.22) 



(4.23) 



(4.24) 



Here yx is the solution of the corresponding spherical field equation for the finite system 
which reads 
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r(¥) 



(47r 



Xd+l)/2 



(d-l)/2 



h 



An 



dx 



+ ^ + — / — {Attx 



.-{d-l)/2 



1 - ( 1 + 2i? ( — 

X 



X 
d~d'' 



exp {-xyx] 



(4.25) 



The expressions ( [4.20| ) and ( |4.21| ) for the infinite system, and ( [4 .241 ) and ( [4.25| ) for the 



finite one, represent actually the verification of the analog of the Privman-Fisher hypothe- 



sigiiJ for the finite-size scaling form of the free energy, formulated for classical systems, for 
the case when the quantum fiuctuations are essential. Note, that in that case one has both 
finite space dimensions and one additional finite dimension that isproportional to the in- 
verse temperature. According to the finite-size scaling hypothesisio one has to expect that 
the temperature dependent scaling field multiplying the universal scaling function will be 
with exponent p = l + rf/z, where the dynamic-critical exponent z expresses the anisotropic 
scaling between space and "temperature" ("imaginary-time") directions. For the system 
considered here z = 1, in full conformity with our results. It seems important enough to 
emphasize that in the quantum low-temperature case there is one nonuniversal pre-factor 
(for the system considered here it is JX) that multiplies the universal finite-size scaling 
function of the free energy, which is not the case for the classical systems (provided their 
free-energy density is normalized per 

Now we pass to the analysis of the spherical field equations for both the finite and infinite 
systems following Ref.0 



V. ANALYSIS OF THE SPHERICAL FIELD EQUATIONS 

A. The infinite system 

After performing the integration in the right-hand side of the equation for the spherical 
field ( [4.21|) and using the integral representation ( |A2|) we get 



1 1 

A ~ a7 



(47r)(''+i)/2 



r 



1 -d 



-l)/2 



+ 



(47r)('i+i)/2 



0(^-l)/2^ 



d-l A 1 
2 ' 2t 



/2 



+ 



02' 

(5.1) 



where 



lC{u,y) = 2Y,{ym)-''K,{2my) 



m=l 



(5.2) 



and Kjj{x) is the other modified Bessel function. 

Eq. ( p.l| ) was derived in Ref. |T by a different technique. The approach used here is 
based on the estimations (^l9|, [4.10|) and has the advantage that it can easily be generalized 
for other boundary conditions.!^ 

In the remainder of this section we will study the effect of the temperature on the 
susceptibility and the equation of state near the quantum critical fixed point. 
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1. Zero-field susceptibility 



After making the field h vanishes, from Eq. ( p.lD we find that the normahzed zero-field 
susceptibility x = on the line A = Ac {t 0~^) is given by 



X 



where yo is the universal solution of 

1 -d 



r 



X2 



2/C 



(5.3) 



(5.4) 



The behavior of the universal constant yo as a function of the dimensionality d of the system 
is shown in FIG. |I]. 

Eq. (|5.3|) tells us that at low-temperature the susceptibility increases as the inverse of 
the square of the temperature above the quantum critical point. 

In what follows we will try to investigate Eq. ( |5.1j ) for different dimensions {1 < d < 3) 
of the system and in different regions of the {t, A) phase diagram. To this end we will 
consider first the function /C(z/, y). For d 7^ 2 its asymptotic form is (for y <^ 1)0 (see also 
Appendix 0) 



d 













d\ 


^ r 1 




2 V 


2) 



-1 



d 



ad-l)y 



l-d _ ^-p 

2 



l-d 



(5.5) 



where ({x) is the Riemann zeta function. Introducing the "shifted" critical value of the 
quantum parameter due to the temperature by 



1 



1 



Ac(t) Xc 27r{'^+i)/2 \X,{t) 



t 



d-l 



r 



d-l 



or 



A?(t) 



1 



Ac ' 27r('^+i)/2 V A. 



d-i 



d- 1 



ad 



\ad-i)\., 



(5.6a) 



(5.6b) 



one has to make a difference between the two cases d < 2 "sign — " and d > 2 "sign +". 
In the first case (1 < < 2), it is possible to define the quantum critical region by the 
inequality 



1 1 
A ~ X 



1 



2vr(<i+i)/2 ^Ac(t) 



t 



d-l 



T 



d-l 



\ad-i)\. 



(5.7) 



For d < 2 the function /C(z/, ?/) 
(outside of the quantum critical region) 



,-1 



by substitution in ( ^.1|) we obtain for A < Ac 
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X 



|r(i-i)| A, 

(47r)"'/2 A, - A 



-1 2/(d-2) 



^2/{d-2)_ 



(5.8) 



We see that the susceptibihty is going to infinity with power law degree when the quantum 
fiuctuations become important (t — > 0"*") and there is no phase transition driven by A in the 
system for dimensions between 1 and 2. 

In the second case (2 < d < 3), after the insertion of Eq. ( pTsD in Eq. ( ^.1|) one has 



X 



|r(i-i)| Xcit) 

(47r)'^/2 A-Ae(t) 



-1 2/{d-2) 



£2/{d~2)^ 



(5.9) 



as a solution for A less than Ac and greater than the critical value Ac(t) of the quantum 
parameter. Here for finite temperatures there is a phase transition driven by the quantum 
parameter A with critical exponent of the (i-dimensional classical spherical model 7 = 
This however is valid only for very close values of A to A(t). For A < Ac(t) the susceptibility 
is infinite. 

In the region where A > Ac the zero-field susceptibility is given by 



X 



(47r: 



,(d+l)/2 



r(i^) VA Ac 



2 



(5.10) 



This result is valid for every d between the lower and the upper quantum critical dimensions. 

Eq. ( [5.5| ) is derived for dimensionalities d ^ 2. The important case d = 2 requires special 
care. In this case Eq. (|5.1| ) takes the form {h = 0) 



A 



1/2 



4:71 



2lT 



1 — exp 



This equation can be solved easily and one gets 



a/2 



t VAc 



For the susceptibility, Eq. ( |5.11|) yields 

A2 



X ~ ^ exp ( 47r 



Ac- A 



tXr 



(5.11) 



(5.12a) 



for 



2lT 



A, 



- 1 



3.1114... 



^ 1 and A < Ac i.e. in the renormalized classical region. For A = Ac 



1 /A 



^ 02 V t 



(5.12b) 



where the universal constant 
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(5.12c) 



was first obtained in the framework o: 
one finite dimension .111 Finally for ^ 
region 



the 3-dimensional classical mean spherical model with 
— 1 3> 1 and A > Ac i.e. in the quantum disordered 



X 



XXr 



_An{\ - Ac 



2tXr 



7r(A - Ac) 



exp 



An 

— (A -Ac 

tXr 



(5.12d) 



The first term of Eq. (|5.12d|) is a particular case of Eq. ( p.lOj ) for d = 2. 

From Eqs. ( |5.12| ) one can transparently see the different behaviors of x{T) in tree regions: 
a) renormalized classical region with exponentially divergence as T — > 0, b) quantum critical 
region with x{T) ~ and crossover lines T ~ |A — Ad, and c) quantum disordered 
region with temperature independent susceptibility (up to exponentially small corrections) 
as T — > 0. The above results ( [5. Ill ) and ( [5. 12] ) coincide in form with those obtained in 
Refs. p|,p!0|for the two dimensional QNLaM in the n ^ oo limit. The only differences are: in 
Eq. (|5.11|) the temperature is scaled by A, and the critical value Ac is given by Ac = 1/Wd(0) 



while for the QNLaM it depends upon the regularization scheme. It will be useful to clarify 
the bulk critical behavior of x found above in the context of FSS theory. First, from Eqs. ( ^.6|) 
one can see that the exponent ip characterizing the shift in the quantum parameter Ac by 
the temperature t is equal to = — 1 in accordance with FSS prediction. Second, let us 
consider the critical behavior of {d+ l)-dimensional classical model with L,- ~ 1/T playing 
the role of a finite-size in the imaginary time direction, i.e. with slab geometry 00°' x L^-. 
FSS calculation (see for example Ref. ^) for the susceptibility at the critical point (in the 
case under consideration this is A = Ac) gives: (i) for 1 < d < 2,0 

x(Ac)~L^^/^\ 

where 7a = 2ux = ^ and x(Ac) = T'^t; 7^ = 2, (ii) for 2 < c/ < 3, 



X(A) ~4^^'^^)/^^|A-A(L, 



-7a 



where the exponent 7a = is pertaining to the infinite size (space) dimensions d. In the 
last case for the shift of Ac we take from Eqs. (|5^ ) Ac — A(L^) ~ f^'^ and again x(Ac) ~ t^'^T, 
So the different types of critical exponent depending on A and T are related by 7^. = 7a/i^a 
or 7j, = ip'-^x {ip = l/uxis the crossover exponent). 



2. Equation of state 

The equation of state of the Hamiltonian (|2.1|) near the quantum critical point is obtained 
after substituting the shifted spherical field (p by the magnetization Ai through the relation 

M = ^, (5.13) 
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in Eq. ( |5.1| ), which allows us to write the equation of state in a scaling form 



1 + 



_6X_ 



(4vr) 



-(d+l)/2 



h 



1/7 



r 



d 



2/C 



d-1 X fM"'^ 
~2~'2 \ t 



(5.14) 



where 



We conclude that near the quantum critical point Eq. ( 5.14|) may be written in general forms 



as 



(5.15a) 



or 



M 



X 



-Ply 



fM{SXM-^^f^,hM-^) 



(5.15b) 



In Eqs. (|5.15| ) fh{x,y) and fM{x,y) are some scaling functions, furthermore 7 = ^^j, 
u = (3 = \ and 5 = are the familiar bulk critical exponents for the {d + 1)- 
dimensional classical spherical model. Eqs. ( [5.15| ) are a direct verification of FSS hypothesis 
in conjunction with classical to quantum critical dimensional crossover. They can be easily 
transformed into the scaling form (Eq. (21)) obtained in Ref. |T^, however here they are 
verified for 1 < d < 3 instead of 2 < (i < 3 (c.f. Ref. p!7|) , i.e. the non-critical case is 
included. 

Hereafter we will try to give an explicit expression of the scaling function fh{x,y) {x = 
SXAi'"^, y = [t/Xy^^Ai"'^) in the neighborhood of the quantum critical fixed point. This 
may be performed, in the case jy^h/Ad ^ 1, with the use of the asymptotic form from 
Eq. (|3D to get (for rf^ 2), 



M J 



(d-2)/2 



1 



2Tr{d+i)/2 \ X 



d-l 



r 



d-i 



C(rf-i) = o. 

(5.16) 



From this equation we find for the scaling function the result 



fh{x,y) 



(47r 



\d/2 



1 + X + 



1 



2T^{d+l)/2- 



r{i-iy 

For the special case d = 2 Eq. (|5.14|) reads 

t 



d 



2/(d-2) 



(5.17) 



6X 



A-kM^I^ + 27rA 



In 



1 — exp 



(5.18) 



which yields 
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h = M^fh{x,y) 
where the scahng function is given by the expression 



fhix,y) = Ay^ 



arcsinh- exp ( 2tt — '^—^ 

2 y 



(5.19) 



(5.20) 



At X = and y ^ 1 (fixed low-temperature and h 0^) the scahng function ( ^.20[ ) 
reduces to 



fhiO,y) ^ y^exp ( . 



(5.21) 



In the region x < — 1, and for y ^ 1 (fixed weak field and t ^ 0+) the corresponding scaling 
function is 



fh{x,y) ^ y^exp An 



x + 1 



(5.22) 



and in region x > — 1 and y ^ 1 we have 



fhix,y) ^ 167r2(x + l)' 



7r(l +x) 



exp —An 



l + x 



(5.23) 



This identifies the zero temperature {y = 0) form of the scaling function ( [5.23| ) with those 
of the 3-dimensional classical spherical model. 



B. System confined to a finite geometry 

When the model Hamiltonian (|2.1|) is confined to the general geometry L'^~'^'xoo'^'xLt-, 
with < d' < d, equation ( [4.25| ) of the spherical field (p takes the form (after some algebra) 



A A, 



-(4vr) 



-{d+l)/2 



r 



1 - d 



-l)/2 



(f){d~l)/2 



l)/2 



{d-l)/2 



0i/2{(Am/t)2 + (L|Z|)2}^/' 



^^''> m,i{d-d') UV2{(Am/t)2 + (L|Z|)2}^/2 



(d^l)/2 ' 02- 



(5.24) 



where 



\l\ = {11 + 11 + ■■■ + ll-d) 



1/2 



and the primed summation indicates that the vector with components m = li = I2 
Id-d' = is excluded. 
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1. Shift of the critical quantum parameter 



The finite-size scaling theory (for a review see Ref. ^) asserts, for the temperature driven 
phase transition, that the phase transition occurring in the system at the thermodynamic 
Umit persists, if the dimension d! of infinite sizes is greater than the lower critical dimen- 
sion of the system. In this case the value of the critical temperature Tc(oo) at which some 
thermodynamic functions exhibit a singularity is shifted to T'c(L) critical temperature for a 
system confined to the general geometry L'^~'^ xoo*^ , when the system is infinite in d' dimen- 
sions and finite in {d — d') -dimensions. In the case when the number of infinite dimensions 
is less than the lower critical dimension, there is no phase transition in the system and the 
singularities of the thermodynamic functions are altered. The critical temperature Tc(cxd) 
in this case is shifted to a pseudocritical temperature, corresponding to the center of the 
rounding of the singularities of the thermodynamic functions, holding in the thermodynamic 
limit. 

In our quantum case, having in mind that we have considered the low-temperature 
behavior of model ( p.l|) in the context of the FSS theory it is convenient to choose the 
quantum parameter A like a critical instead of the temperature t and to consider our system 
confined to the geometry L'^~'^xx)'^xL^. So the shifted critical quantum parameter Ac(t, L) = 
XtL is obtained by setting = in Eq. (p.24|). This gives 



1 



1 



r(¥) 



A, 47r('^+i)/2 



j2 {{XtLm/tr + {L\i\r) 



2^{l-'^)/2 



(5.25) 



m,l{d-d') 

The sum in the r.h.s. of Eq. (|5.25| ) is convergent for d' > 2, however it can be expressed 



in terms of the Epstein zeta function (see Ref. ^5]) 



Z 




l-d 
2 



m 



(5.26) 



which can be regarded as the generalized [d — d' + l)-dimensional analog of the Riemann 



zeta function C(-^)- In the case under consideration the Epstein zeta function has only a 
simple pole at d' = 2 and may be analytically continued for < < 2 to give a meaning to 
Eq. ( |5.25|) for rf' < 2 as well. It is hard to investigate the sum appearing in Eq. ( |5.26| ). The 



anisotropy of the sum L'^lf+- ■ + ( j) is an additional problem. For these reasons 

we will try to solve it asymptotically, considering different regimes of the temperature, 
depending on that whether L -C ^ or L ^ ^ which will be called, respectively, very 
low-temperature regime and low-temperature regime. To do that it is more convenient to 
write Eq. ( |5.25| ) in the following form 



11 1 v-V ±± I r2i2 / ^tL\ 2 



i,l{d-d') 

This will be the starting equation for the two cases we will consider below. 



(5.27) 
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a. Low-temperature regime ^ <^ L In this case we single out, from the sum in 



Eq. ( |5.27|) , the term with Z = to obtain 

d-l 



1 



1 



E' E 

l{d-d')m=-co 



Cid-l] 



dxx 2 exp 



-X 



t 



m 



(5.28) 



After the apphcation of the Jacobi identity ( |4.6|) to the sum over m and the calculation of 
the arising integrals in the resulting expression we obtain the final result 



1 



1 



1 



t 



27,(d+l)/2 \XtL ^ 



d-l 



d- 1 



C(rf-i) + 



t L 



2-d 



XtL 47r'^/2 



r 1-1 EV 

^ ^ l{d-d' 



-d 



+ 



TT 



E E 

i(d-d') m=l 



(i/2-1 



t 



2n- — Lm\l\ 



(5.29) 



The first term of the r.h.s of Eq. ( p.29| ) is the shift of the critical quantum parameter 
(see Eq. ([5. 6])) due to the presence of the quantum effects in the system. The second term 
is a correction resulting from the finite sizes. It is just the shift due the finite-size effects in 
the (i-dimensional spherical mo delS multiplied by the temperature scaled to the quantum 
parameter. Here the {d — (i')-fold sum may be continued analytically beyond its domain 
of convergence with respect to "c?" and "d"' (which is 2 < d' < d). The last term is 
exponentially small in the considered limit i.e. ^ -C L. 

In the borderline case d = 2, one can see that in the first and second terms in the r.h.s 
of Eq. ( |5.29| ) singularities take place. In Appendix ^ we can show that they cancel and 
Eq. (|]29D then yields 



1 



1 

a; 



2nX 



tL 



5o + 7ij + In ^ + 2 f; i^o (2n^Lm\l\ \ [ , (5. 

*^ l{2-d')m=l \ tL /J 



30) 



where 7^ = .577... is the Euler constant. 

In the particular case of strip geometry d' = 1, the shift is given by (see Appendix 



1 



1 

a: 



27rA 



tL 



7i? + ln 



tL 



47rA 



tL 



^^EEA-oK-')}- 



(5.31) 



For the fully finite geometry case {d' = 0) we obtain for the shift (see Appendix 



1 



t 



Ac 27rA 



tL 



, tL 



In MI + 2 V' E ^0 { 2n^LmJq7q] 

^ V A,, V 1 V 



(5.32) 



Let us note that the last term in the r.h.s of Eqs. ( |5.31| ) and ( ^.32|) gives exponentially small 
corrections and in the considered limit they can be omitted. 
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b. Very low-temperature regime j ^ L To get an appropriate expression now, we single 
out the term corresponding to m = from the sum in Eq. ( |5.27| ). This leads to 



1 



1 



l{d~d') 



-d 



l{d-d') m=l 



d-3 

ax X 2 exp 



-X 



lH' + ( --^ 



X 



tL 



m 



(5.33) 



The next step consists in the application of the Jacobi identity to the {d — c?') -dimensional 
sum and the calculation of the arising integrals to obtain 



A 



tL 



1 



A. 47r('^+i)/2 



l{d-d') 



+ 



^l/2+d72-d / ^ 



TT 



2vrK+i)/2- 



d'-i 



d'-i 



i- ] ad' - 11 



l(d-d') m=l 



(5.34) 



Here, in the r.h.s, the first term is the expression of the shift of the critical quantum pa- 
rameter, at zero temperature ,0 due to the finite sizes of the system. This is equivalent to 
the shift of a (c? + l)-dimensional spherical model confined to the geometry xoo*^ . 

The second term gives corrections due to the quantum effects. This is the shift of criti- 
cal quantum parameter of a c/'-dimensional infinite system multiplied by the volume of a 
{d — c?') -dimensional hypercube. The third term is exponentially small in the limit of very 
low temperatures. The singularity for rf' = 1 of the r.h.s of Eq. ( ^.34| ) is fictitious. So for 
d' = 1 Eq. (|5.34|) yields (see Appendix 



1 1 

XtL Ac 



l~d 



2n 



In 



l{d-l)m=l 




(5.35) 



where 



In the particular case d' = 1 and d = 2 the constant Cq is given by Cq = — In An and 

(5.36) 



1 



1 

a: 



II ^1 

7e + In ■ 



2nL 



A-KtL 



+4f:f: A„(2.^™A) 

1=1 m=l ^ ^ J 



Comparing between Eqs. ( ^.31 ) and ( 5.36| ) one can see the crucial role (in symmetric 
form) of L or Xti/t in the low-temperature regime and the very low-temperature one, 
respectively. 

In another particular case of a two-dimensional bloc geometry d' = and d = 2, using 
the Hardy formula (see Appendix P) one gets 
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(5.37) 

Instead of the previous case of low-temperature regime here the lower quantum critical 
dimension d' = 1 is responsible for the logarithmic dependence in Eq. ( [5.351 ). This is the 
reason for the significant difference between Eqs. ( |5.36|) and ( ^.37|) . 

The above obtained equations for Xti will be exploited later for the study of the two 
dimensional case. 



2. Zero-field susceptibility 
It is possible to transform Eq. ( |5.24|) in the following equivalent forms 



(L0 



l/2\d-l 



(47r)('^+i)/2 



d-d' + 1 



where 



I + a^m^ 



or 



l-d 



5\ + 



X \t 



(A0VVt)d-i 

—2K,tL ( 

- V 2 



r 



1 - d 



d-d' + 1 



X(p 



1/2 



2t 



where 



ICa{jy\p,y) = ICi{iy\p;ay) = — xi 7^, y > 0. 



(5.38a) 



l/>0, V = li + li + --- + l' (5.38b) 



(5.39a) 



(5.39b) 



The functions /Ca(z/|j9;y) and ACa(z/|p;y) are anisotropic generalizations of the /C-function 
introduced in Ref. 

Eqs. ( 5.38 ) and ( ^.39 ) show that the correlation length ^ = 0"^/^ -^wi scale like 



(5.40a) 



or like 



19 



'T'Ha 



(5.40b) 



which suggests also that there will be some kind of interplay (competition) between the 
finite-size and the quantum effects. 

Hereafter we will try to find the behavior of the susceptibility x = 0"^ as a function 
of the temperature t and the size L of the system. For simplicity, in the remainder of this 
section, we will investigate the free field case {h = 0). 

1. For <C 1, after using the asymptotic form of the function defined in ( |5.38b|) (see 

Appendix 0) 



7]-P/2 /„ 



1. 



^,-C7„(p|.)--r(-.), 



(5.41) 



Eq. ( p:24D reads (d^ 2, 1 < < 3) 

{d'-2)/2 



J. Td'-d 



A (47r)'^'/2 



d' 
2" 



+ 



1 



4^{d+l)/2 



r 



d-l 



m,l{d-d% 



^m) +{Llf 



l-d 
2 



= 0. 

(5.42) 



Now we will examine Eq. ( |5.42|) in different regimes of t and L and for different geometries 
of the lattice: 

l.a. jcj)^^"^ <^ 1 and ^ ^ 1: In this case Eq. ( ^.421) transforms into (up to exponentially 
small corrections c.f. with Eq. ( [5.29 )) 



„ t L^ -^ ^ / d' 

= SX + -- — -^r 1 

A {AtiY'/^ V 2 



XiTt'ft I 2 



r I' 

^ l(d-d') 



d'-2)/2 



2-d 



+ 



1 



d~l 



r(^U(rf 



(5.43) 



This equation has different types of solutions depending on whether the dimensionality d is 
above or below the classical critical dimension 2. 

At A = Ac and when d' < 2 < d < 3 (i.e. when there is no phase transition in the system) 
we obtain for the zero-field susceptibility 



X 



^ \ -2 /^j^\2{d-d')/{2-d') 



Ac 



A. 



r(¥) 



7,(d-d'+i)/2 r (1 - f ) 



2 

2-d' 



(5.44) 



However, for 1 < d < 2, Eq. ( 5.43 ) has no solution at A = Ac obeying the initial condition 

Eq. ( ^.44| ) generalizes the bulk result ( |5.3| ) for d close to the upper quantum critical 
dimension i.e. d = 3. 

At the shifted critical quantum parameter Ac(t) given by Eq. ( |5.6| ) we get 
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X 



2^'-2 r(f-i)_ , 



2-d 



^{d-d')/2 r (1 - f ) ^ 



(d-d') 



2 

2-ii' 



(5.45) 



However this solution is valid only for 3 > (i > 2 > i.e. here again there is no phase 
transition in the system. 

l.b. < 1 and f < 1: In this case, Eq. ( ^.421) gives (up to exponentially small 

corrections c.f. with Eq. (|5.34|) ) 



J- Td'-d 

6X + --^T 



1-^ 

2 



,K-2)/2 



l-d 



+ 



4j^{d+l)/2 



d- 1 



l{d-d') 



'd'-d 



r 



d' -I 



d'-i 



C(rf'-i). 



(5.46) 



Here we find that the solutions of Eq. ( ^.46|) depend upon whether the dimensionality d' <1 
or d' > 1. 

At A = Ac and for 1 < d' < 2, Eq. (|^ has 



X 



tL J 



2/(2-d') 



(d-d'+l)/2 Y (-\ _ d[\ I*"! 

^ 2)l{d-d') 



2 

2-d' 



(5.47) 



as a solution. For < d' < 1, however, it has no solution obeying the initially imposed 
restriction ^0^/^ < 1. 

At the shifted critical quantum parameter Xc{L) given by0 



1 
A 



-i-d 



r 



d-i 



l{d-d') 



-d 



Eq. ( |5.46|) has a solution obeying the initial condition 
this case the susceptibility behaves like 



X 



_J_K 



2 r 



^ 1 only for (i' = 1 + £ and in 



(5.48) 



2. For L</)^/^ ^ 1 from Eqs. ( 5.38 ) and Eq. ( |A9| ) we get once again Eq. ( |5.43| ). In spite 



of the fact that we have the same equation as in the case 



<^ 1, the expected solutions 



for the susceptibility may be different because of the new imposed condition. Here also we 
will consider the two limiting cases of low-temperature and very low-temperature regimes. 

2. a. Lcf?-/'^ <^ 1 and ^ ^ 1: In this case Eq. ( ^.42| ) again is transformed into Eq. (|5.43|) 
and we obtain at A = Ac the solution given by Eq. ( |5.44| ), which is valid only for d' < 2 < 
d < 3, i.e. we have the same solution as in the previous case i.e. j^^''^ ^ 1. 

At A = Xc{t), we formally obtain Eq. ( p.45| ) which, however, may be considered like a 
solution only in the neighborhood of the lower classical critical dimension d = 2. For the 
cylindric geometry {d' = 1 and ci = 2 + e) we get 
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X 



1- 2(^i^-1^47r) 



(5.49) 



This result is contained in Eq. (30.109) of Ref. ^ in the large n-limit case for the NLQcrM. 
In the case of slab geometry d — d! = l{d = 2 + e,d' = l + e) instead of ( ^.49|) we obtain 



X 



[Tie] 



[l-£(7^-ln2) -£ln£] 



(5.50) 



In the case of a bloc geometry {d = 2 + e and d' = 0) we find the following behavior for 
the susceptibility 



X 



2ne 



1-4 |7i.-ln 



Mil 

471-2 



(5.51) 



For the case of "quasi-bloc geometry" {d = 2 + e and d' = e) we get 



X 



2Tie 



ne 



- 27^ + ln— -21n 



Mil 

2v^ 



(5.52) 



Here the appearance of e in the denominator in formulas ( |5.48| )-( |5.52| ) signalizes that the 
scaling in its simple form will fail at e = 0. 

2.b. L(j)^/'^ < 1 and x ^ 1- Here we find that Eq. (|5.46| ) is valid, and it has Eq. (^.47|) 
as a solution at A = Ac and for < rf' < 1. For 1 < rf' < 2 the susceptibility is given by 



X 



2 r(^) 



V2r(i-f) 



7r 



U{d' - 1) 



2 

2-d' 



(5.53) 



At the shifted critical point \c{L), for the susceptibility we obtain Eq. (|5.53| ) under the 
restriction 2 > c?' > 1, which guarantees the positiveness of the quantity under brackets. 

When A < Ac for 1 < ci < 3 and d' < 2, i.e. when there is no phase transition in the 
system, we obtain 



X 



(47r 



r(i-f) 



2 

2-d' 



^-2/{2-d') j2[d-d')/(2-d') 



(5.54) 



If > 2 there is a phase transition in the system at the shifted value of the critical 
quantum parameter XtL (the shift in this case is due to the quantum and finite-size effects) 
and Eq. ( p.42| ) transforms to 

A 



1 - 



\tL ^(47r)^/2 
which has the following solutions 



L'^'-'^ ( d' 
t- -TT^rr 1 



M'-2)/2 



X 



(47r)''V2 



1 



A 



2_ 

2-d' 



00, 



^-2/(2-d')]^2{d-d')/{2-d')^ A > AtL 

A < \tL 



(5.55) 



(5.56) 



Let us notice that Eqs. ( ^.54| ) and ( ^.56|) are the finite-size forms, for the susceptibility, of 
Eqs. (|5.^ ) and ( |5.9| ), respectively, found for the bulk system. 
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3. Two-dimensional case 



The two-dimensional case needs special treatment because of its physical reasonability 
and the increasing interest in the context of the quantum critical phenomena.i'i^El From 
Eq. ( |5.24|) for d = 2 and in the absence of a magnetic field h = we get 



m,l{2-d') 



1/2- 



01/2 1 . 

17 47r ^ . (^^2 + ^2^2)1/2 



(5.57) 



Introducing the scaling functions Yf = ^0^/^ and Yf = Ltp^/"^^ where the superscript d! 
denotes the number of infinite dimensions in the system, and the scaling variable a = ^ it 
is easy to write Eq. ( p.57[ ) in the scaling forms given in Eqs. ( p.38| ) and ( |5.39| ). The solutions 



of the obtained scaling equations will depend on the number of the infinite dimensions in 
the system. Here we will consider the two most important particular cases: strip geometry 
d' = 1 and bloc geometry d' = 0. Our analysis will be confined to the study of the behavior 
of the scaling functions at the critical value of the quantum parameter Ac, and at the shifted 
critical quantum parameter (see Section |V B 1| ) . It is difficult to solve Eq (|5.57|) by using 



an analytic approach that is for what we will give a numerical treatment of the problem. It 
is, however, possible to consider the two limits: a ^ 1 i.e. the low-temperature regime and 
a -C 1 i.e. the very low-temperature regime. 

Strip geometry (d' = In this case in the r.h.s of Eq. ( |5.57|) we have a two-fold sum 



which permits a numerical analysis of the geometry under consideration. FIG. graphs the 
variation of the scaling functions Y^ and Yl against the variable a at A = Ac. This shows 
that for a comparatively small value of the scaling variable a ~ 5 the finite-size behavior 
(see the curve of the function Y^{a)) merges in the low temperature bulk one, while the 
behavior of Yl{a) shows that for relatively not very low-temperatures (a ~ |, L- fixed) the 
system simulates the behavior of a three-dimensional classical spherical model with one 
finite dimensions. The mathematical reasons for this are the exponentially small values of 
the corrections, as we will show below. 

Bloc geometry (d' = Oj: In this case the three-fold sum in the r.h.s of Eq. ( |5.57| ) is 



not an obstacle to analyze it numerically. For A = Ac the behavior of the scaling functions 
Y^(a) and Y^(a) is presented in FIG. 0. They have the same qualitative behavior as in the 
strip geometry the only difference is the appearance of a new universal number for t = 0, 
i.e. f2, instead of the constant 6 as a consequence of the asymmetry of the sum in the 
low-temperature and the very low-temperature regimes. 

Analytically for arbitrary values of the number of infinite dimensions d' we can treat 
first the problem in the low-temperature regime (a ^ 1). In this limit Eq. (|5.57|) can be 
transformed into (up to small corrections (9(e~^'^")) 

5\ = ^ln2sinh-0i/2 _ J_ y' ^ (L(j)'/^\l\) . (5.58) 
27rA 2t^ 27rA ^ ^ ^ ' ^ ^ 

l{2~d') 

For A = Ac Eq. (|5^) has the solut ion 
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,-V.^J_e^(2_,0;|(^)'^.p(-fe) (5.59, 

i.e. the finite-size corrections to the bulk behavior are exponentially small. 

In the very low-temperature regime (a <^ 1), Eq. ( |5.57| ) reads (up to 0{e^'^'^^"')) 

^^ = 17-172^ \l\ + 7^^2-J l^K^[j4>^ , (5.60) 

l{2-d') ' ' ^ ^ m=l ^ / 



which has the solutions: 



1 ^ /27r f LXrV^ ( Xr 



for d' = 1, and 



-"4"4{2k4i:'("= + ^'^'T^4 -p(-4) (5.62) 



for = 0. 

In Section |VB 1| an analytic continuation of the shift of the critical quantum parameter 
for (i = 2 was presented. It is possible to consider the solutions of Eq. ( |5.57| ) at A = XtL 
(from Eqs. (|5.31|) , ( |5.32|) , ( ^.36| ) and ( ^.371) ) and for different geometries. In this case the 
scaling functions F/, Yl, and Y^ are graphed in FIG. ^. For d' = 1 again we see that a 
symmetry between the two limits a -C 1 and a ^ 1 take place, since the scaling functions 

and Yl are limited by the universal constant S. The asymmetric case d' = has two 
different constants and T^l, limiting the solutions of Y^ and Yl from above. 

The constants S, and T^l are obtained from the asymptotic analysis (with respect to 
a) of Eq. ( p37D for X = XtL- 

In the limit a ^ 1 for arbitrary values of d' we get (from Eq. (|5.58|) ) 

i?o + 7i^ + ln^= 5^Vo(L0i/2|Z|), (5.63) 

l{2-d') 

where the equation of XtL from Eq. ( ^.3Q| ) is used. Eq. ( |5.63|) has the solutions: 

'--"'-{It i'Zl: («-^*) 

where the universal numbers S = 7.061132... and T^l = 4.317795... are the solutions of the 
scaling equation ( |5.63|) for d' = 1 and d' = 0, respectively. 

In the opposite limit a <^ 1, for d' = 1, we get from Eqs. ( |5.36| ) and ( p.6U| ) the equation 

7e + In = 2 > : i^o ( ) , (5.65) 



Ant 



2 Y: (^<l>'^'rn 

m=l ^ 
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which has 



A 



!L -1/2 

t ^ 



universal solution. For rf' = we have 



A 



exp 



t 



A 



'-^6"^ -6 = 



(5.66) 



(5.67) 



obtained from Eqs. ( |5.37| ) and ( [5.60| ), where we have used the identity ([All| ) 
From Eq. (|5.67|) we obtain the universal result 



A 



-1/2 
t ^ 



St = 6.028966... 



(5.68) 



We finally conclude that if we take A = Ac the scaling functions Yf and Yf have a similar 
qualitative behavior weakly depending on the geometry (i.e. bloc rf' = or strip d! = 1) 
of the system. However, for a given geometry one distinguishes quite different quantitative 
behavior of the scaling functions depending on whether the quantum parameter A is fixed 
at its critical value, i.e. A = Ac, or takes "running" values XtL obtained from the "shift 



equations" (OTl), (OD, (KM or (K37\) 



4- Equation of state 



The equation of state of the Hamiltonian (|2.1|) for dimensionalities 1 < d < 3 is given by 
(see Eqs. ( ^AM and j^M)) 



= 6\- (47r 



,-(d+l)/2 



r 



1 -d 



1^ \ (d-l)/2 

M 



{d-l)/2 



(5.69) 



It is straightforward to write this equation in a similar form as in Eq. ( [5.38| ) or Eq. ( p. 391 ), 
h = M'ff: ^6\M-^'^ y , L-^/'^M-^/^^ , (5.70a) 



I.e. 



or 



(5.70b) 



Eqs. ( |5.7U| ) are generalizations of Eqs. ( 5.15 ) in the case of systems confined to a finite 
geometry. The appearance of an additional variable ^ is a consequence of the fact that the 
system under consideration may be regarded as an "hyperparallelepiped" (in not necessary 
an Euclidean space) of linear size L in d — d' directions and of linear size in one direction 
with periodic boundary conditions. 
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VI. SUMMARY 



Since exact solvability is a rare event in statistical physics, the model under consider- 
ation yields a conspicuous possibility to investigate the interplay of quantum and classical 
fluctuations as a function of the dimensionality and the geometry of the system in an exact 
manner. Its relation with the QNLaM in the n ^ oo limit may serve as an illustration of 
Stanley's arguments of the relevance of the spherical approximations in the quantum case. 
Let us note, however, that the use of such types of arguments needs an additional more 
subtle treatment in the finite size case. For this reason we gave a brief overview of the bulk 
low-temperature properties, which are similar to those obtained by saddle-point calculation 
for the QNLaM (see Section [VAl) . 

The Hamiltonian ( ^.11 ) can be obtained also in the "hard-coupling limit" from a more 
realistic model that takes into account the quartic self-interaction term Qj, by the ansatz 
Ql ^ J/ J2e Qi frequently used in the theory of structural phase transitions (see Refs. 
and Section |T|). 

An analog of the Privman-Fisher hypothesisiS for the FSS form of the free energy in 
the quantum case was shown to be consistent with the exact results obtained in Section [rV| . 
Furthermore, the explicit expression for the universal scaling function is presented and as a 
consequence the scaling form of the equation for the spherical field. 

So the Hamiltonian ( p.l|) can be thought of as a simple but rather general model to 
test some analytical and numerical techniques in the theory of magnetic and structural 
phase transitions. The discussion of the results obtained, in Section |V A| , serves as a basis 
for further FSS investigations. Identifying the temperature, which governs the crossover 
between the classical and the quantum fluctuations as an additional temporal dimension 
one makes possible the use of the methods of FSS theory in a very effective way. 

In Subsection |V B 1| the shift of the critical quantum parameter A as a consequence of the 
quantum and finite-size effects is obtained. In comparison with the classical case (for details 
see Ref. and references therein) here the problem is rather complicated by the presence of 
two finite characteristic lengths L and Lr- We observe a competition between finite size and 
quantum effects which reflects the appearance of two regimes: low-temperature {L ^ L^) 
and very low-temperature {L ^ L^-). The behavior of the shift is analyzed in some actual 
cases of concrete geometries e.g. strip and bloc. 

In the parameter space (temperature t and quantum parameter A), where quantum 
zero point fluctuations are relevant, there are three distinct regions named "renormalized 
classical", "quantum critical", and "quantum disordered". The existence of these regions 
in conjunction with both regimes: low-temperature and very low-temperature make the 
model a useful tool for the exploration of the qualitative behavior of an important class of 
systems. 

In Subsection |V B 2| the susceptibility (or the correlation length) is calculated and the 
critical behavior of the system in different regimes and geometries is analyzed. We have 
studied the model ( ^.11 ) via e-expansion in order to illustrate the effects of the dimensionality 
d on the existence and properties of the ordered phase. An indicative example is given by 
Eqs. (|5.49|) and (|5.50|) , while the former is known (see Ref. ^) the last one is quite different 
and new. These shows that one must be accurate in taking the limit e — > O"*" in the context 
of the FSS theory. 
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In Subsection |VB 3| , special attention is paid to the two-dimensional case. The two 
important cases of strip and bloc geometries are considered. The universal constant G given 
by Eq. ( p.l2c| ), which characterizes the bulk system, is changed to a set of new universal 
constants: fl (see Eq. ( ^.62| )), S and (see Eq. (|5.64|) ), and (see Eq. ( ^.68|) ). The 
appearance of new universal constants reflects the new situation, when there are two relevant 
values of the quantum parameter A: A = Ac in the bulk case and A = XtL in the case of finite 
geometries. Due to their universality these constants may play an important role even in 
studying more complicated models. The behaviors of the scaling functions at the bulk critical 
quantum parameter Ac and the shifted critical quantum parameter Xf^ are given in FIGs. |^ 
and 1^. Lt- is the main characteristic length and the corrections are exponentially small 
in the case of low-temperature regime, and vice versa in the case of very low-temperature 
regime. 

The equation of state, for the system confined to the general geometry L'^~'^ xoo'^ xLr, 
is obtained in Subsection V B4 . This reflects the modifications of the scaling functions as a 
consequence of the finite sizes and the temperature. 

Finally let us note that this treatment is not restricted to the Hamiltonian ( |2.1| ), but 
it can be applied to a wide class of finite lattice models (e.g. directly to the anharmonic 
crystal model see Refs. p3H36|) and it can also provide a methodology for seeking different 
quantum finite-size effects in such systems. 
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APPENDIX A: ASYMPTOTICS OF THE FUNCTIONS /Ca(i^|p,y) 

In this appendix we will sketch a way to find the asymptotic behavior of the functions 
]Ca{iy\p,y) defined in section |V B| (see Eq. ( p.38| )). They have the following form 



lCa{p\p,y)= =y7^, y>0, (Ala) 

m,i(p-i) [yVl +a?rTn?\ 

where 

f = ll + ll + --- + lU. (Alb) 
By the use of the integral representation of the modified Bessel function 

(^2Vit) = K.^ (^2v^) = i (^^y^^ ^ x-''-^e-'''-'/''dx (A2) 
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m,l{p—l) 



and the Jacobi identity for a p-dimensional lattice sum 

m,l(p-^) 

we may write (|A1|) as 

MHP,y) = ^r(|-.)r' 



^(l^+rrP/a'^)/t 



(A3) 



2a 



-y 



-a;(«^+m2/a2) 



'7r\P/2 

a I -1 



(A4) 

\ .7; / 

m,l{p-l) 

Let us notice that the two terms in the square brackets in the last equality cannot be 
integrated separately, since they diverge. Nevertheless, in order to encounter this divergence, 
we can transform further ( |A4| ) by adding and subtracting the unity from exp(— xj/^/tt^), 
which enables us to write down (after some algebra) the result 



jr2u-p/2 1 



TT 



2u- 



p/' r (I - z.) 



2a 



,2u 



m,i(p— 1) L 



9 2 \ ^~pI2 

a? TT^ 



2\ '^-p/2 



(A5) 



where 



Ca(pk) = lim / dxx^P^"^^ 



m,l{j}—l) 



"K\P/2 

al -) 



(A6a) 




2 - I., 5 ( /2 + 



2 \ y-p/2 



dmdP-H- 



1^,5 I' + ^ 



(1' + ^) 



v-p/2 



(A6b) 



is the Madelung-type constant and r[a,a;] is the incomplete gamma function. 

We see from Eq. ( [AS] ) that the shift of the critical quantum parameter is given by the 
Madelung type constant ( |A6|) instead of the sum in Eq. (|5.25|) . Indeed it is possible to 
show that these two representations are equivalent. This may be done, following Ref. by 
starting from the Jacobi identity Eq. (|A3|) , where we multiply the two sides by ^p/^-^-i and 
integrating over 5 to obtain the key equation 



CM") 



E 

m,i(p-l) 



, r 



SPl 



2-v 



p/2-v 



r 



+a7r- 



p/2~2u 



7r2(j2+a2m2y 



I o 



IT 



p/2 



(A7) 
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Finally from Eq. (|A7| ) we easily see easily that the integration constant Ca(p|i^) may be 
written in two different forms. In the first case we take the limit 5 — oo and obtain 



(A8) 



In the other case we take the limit 6^0, and then both first and last terms in the r.h.s of 
Eq. (^) yields Eq. (^). 

Using a similar procedure we find, for the functions ICa{i^\p,y) defined in ( |5.39b|) , the 
following expression 



p/2 



^2u-p/2 r (f - Z/) 



2aP 



-1 



y 



2v 



I' 



m 



m,l{p—l) 

Here the Madelung type constant is given by 



u-p/2 



— + m' 



v-p/2- 



■ (A9) 



CahW) = lim / dxx'^P-"-^ 



m,l{p—l) 



x) 



lim < 

<5^0 I ^ 



in,l{p—l) 



u-p/2 



— oo J — oo 



dmdP-H- 



r 



- 1^,5 [K + m 



iy-p/2 




m,l{p-l) 



(AlOa) 



(AlOb) 
(AlOc) 



Eqs. (|A5|) and (|A9D are slight generalizations (for the anisotropic case a 7^ 1) of the result 
obtained in Ref. H2| from one side, and get in touch with the Watson type sums proposed 



earlier in Ref. 19 from the other (see also Ref. El 



If we put in Eqs. ([A5|) or ([A9|) d = 2, d' = and a = 1 we obtain the identity 



E 



,exp f-yy/ll + lA /1\ ^ /f 



+2nJ2' I 



1 



(All) 
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APPENDIX B: SHIFT OF THE CRITICAL QUANTUM PARAMETER FOR 

SOME PARTICULAR GEOMETRIES 

Our task in this appendix is to explain how the shift of the critical quantum parameter 
in the special cases d = 2 (for the low-temperature regime) Eq. (|5.3CI| ) and d' = 1 (for the 
very low-temperature regime) Eq. ( ^.35| ) are obtained. 



1. Low temperature regime, d = 2 



The expression for the shift is given by Eq. ( |5.29| ) for an arbitrary dimensionality of the 
system. To calculate the shift for c? = 2 we proceed as follows: we first put d = 2 + e and 
expand the obtained expressions around e = 0. So, for the functions in the first term of 
Eq. (|]29D we obtain 

(l-|lnvr)+0(52), 

l+e 



^-(3+e)/2 ^ ^-3/2 



=v^[l-|(21n2 + 7^) +0{e'^ 



where 7^ is Euler's constant. For the second term we obtain 

'e\ 2 



r 



and 



lim y"'\i?-'^ = z 



l{d-d') 



e 



1 - ^lE^ 



+ Oie) 



(Bl) 



Finally after the substitution of the above expansions in Eq. ( [5.29| ) and then taking the limit 

e ^ we obtain Eq. (|CT| ). 

Let us now consider the particular cases d' = 1 and d' = 

The case d' = 1: Here the {d — (i')-dimensional sum in Eq. (|BlD reduces to0 



Y^r' = 2C{e) = -1 - £ln(27r) + 0(e^ 



(B2) 



hence the constant Bq = —\n 2ti and the result ( |5.31| ) emerges. 

The case d' = 0: Now the sum in Eq. ( pi|) is two-dimensional. Then we obtain the 
following result due to HardyEl 
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E' {ll + lir^" = ^ae/2)P{e/2) 



where 



'-1 



^-(2/ + l)^' 

To obtain the expansion of the functions in Eq. (p^ ) we use ( P2|) and 

1 2 



/3(£) = 4-^<'C(^,7)-C(^,7 



4 



2 2 27rv^ ^ ^ 



obtained with the aid of the expansionil 



Inr(a) - -ln27r 
^ ' 2 



for the incomplete zeta function ({s, a). Finally we get 



-1 -eln 



(B3) 



(B4) 



(B5) 



which makes possible the calculation of the shift, Eq. ( p.32| ), of the critical quantum pa- 
rameter for the fully finite geometry. 



2. Very low temperature regime, d' = 1 



In this limit the shift is given by Eq. ( p. 34] ). Here we will set d' = 1 + e and expand 
around e = 0. In the first term of Eq. (|5.34|) we have a (c? — c?') -dimensional sum, which 
can be replaced by the corresponding Epstein zeta function. Then we can use the following 



asymptotic form! 



by 




-d 



z 



l(d~d') 



{l\d-l] 



r(¥) 



Co + 0{e) 



d! = l + e, 



(B6) 



where the expressions for the constants Co are quite complicated except for some special 
cases (see Refs. e.g. for ci = 2, c?' = 1 — £ we have Co = 7^; — ln(47r) (c.f. Ref. Eq. 
(30.104)). 

Further for the function in the second term of Eq. (|5.34|) we use the expansions given in 
the preceding subsection. After substituting the obtained expansions in the basic expres- 
sion (Eq.( p.3^ )) for the shift, we obtain Eq. ( p.35| ). 
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FIGURES 




FIG. 1. The dependence of the universal constant yo upon the dimensionaUty d. The constant 
= 0.962424... is obtained for the two-dimensional system (see Eq. ( |5.12cD ) 
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FIG. 2. The effects of the finite-size geometry on the bulk behavior of cf?'/'^ for the two dimen- 
sional case at A = Ac- The superscript d! in Yf = Lcf)^^"^ and Y^^' = ^<p^^'^ indicates the number 
of infinite dimensions in the system. The scaling variable a = j^. The universal numbers are 
e = 0.962424... (see Eq. (|l2^) and n = 1.511955... 
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FIG. 3. The same as in FIG. ^ but for A = XtL and a = j^. The universal numbers are: 
H = 7.061132..., St = 6.028966... and = 4.317795.... 
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